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Abstract 

We construct the exact time dependent density profile for a superposition of the ground and 
singly excited states of a harmonically trapped one dimensional Bose-Einstein condensate in 
the limit of strongly interacting particles, the Tonks- Girardeau gas. Results of an off resonant 
light scattering experiment probing the system could allow to determine the number of particles 
contained in the g well as the coefficients of the superposition. 
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Density fluctuations in dilute Bose-Einstein condensates were traditionally described by 
using a mean- field theory [l, 2], and a similar approach has also been proposed specifically 
for the explanation of the properties of elongated pencil shaped samples [3( , observed also 
in experiments {4, s|. Theoretical descriptions of such quasi one-dimensional systems 
have used the hydrodynamic approximation and have described the properties of 

the fluctuations as corrections to an approximate static density. Problems, however with 
time dependent mean field theories have been pointed out in [8]. More recent experiments 
reported the confinement of Rb atoms in a quantum wire geometry [s[ [loj], where the 
ratio of the interaction to kinetic energy significantly exceeds unity thus approaching 
the Tonks-Girardeau (TG) limit. These works have turned the purely mathematical 



model considered in classic papers 
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applications. For recent reviews see 14 
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15|- 



13( into a real physical system with potential 
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n the case of strongly interacting bosons, when the TG model can be applied 16, 



one can start from the many body wave function of the system and consider 



the exact time dependence of the density determined by the trapping frequency. These 
space and time dependent oscillations give rise to a modulation of a weak probing field 
which can be observed in principle, and can provide information on the properties of the 
condensate without destroying it. 

We give first a microscopic explanation of the observed oscillations based on the the- 
ory of a strongly interacting, one dimensional, harmonically trapped Bose gas. In this 
framework the ground state wave function of the system is equivalent to a noninteracting 
one dimensional Fermi system [121 ] . Its density oscillations are the consequence of the 
quantum mechanical superposition of the ground state with energy E , and one or a few 
excited states. Considering the simplest case of having only a single excitation with energy 
ftu, where u is the angular frequency determined by the harmonic trapping potential, the 
wave function of the condensate is 



ttsfci, ...,x N ,t) = cbttoe-*^* + c 1 tt 1 e-*<3 l+w > t 



(1) 



where \l/o is the ground state and \&i is the first excited stationary state, both being 
symmetric functions of the particle coordinates: xi, . . . ,xn- This wave function yields the 
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time dependent particle density 



g(x,t) = Nf\fy s (x,x 2 . . . ,x Nj t)\ 2 dx2 ...dx 



N 



(2) 



which will obviously exhibit oscillations with frequency to. The explicit form of these are 



g(x,t) = \co\ 2 Q (x) + |ci| 2 £i(x) + 2RecJco£>oi(^)e 



-iujt 



(3) 



where go and Q\ are the particle densities of the ground and excited states respectively 
while £>oi ( x ) — N J^/q(x,X2 ■ ■ ■ ,xn)^i(x,x 2 ■ ■ ■ ,x^)dx 2 ■ ■ ■ dx^ determines the spatial 
pattern of the density oscillations. 

In order to calculate explicitly the terms in this sum we recall the construction given 

for 
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20] . The ground state N particle wave function can be written as 



(N—1,N) 

det ip n {xj) 



]1 sign(x fc 



l<j<k<N 



where 



(4) 



(5) 



2 n nU^/7T 

is the nth normalized harmonic oscillator eigenfunction, where i = (h/mcu) 1 ^ 2 , m is the 
mass of a particle, while the product of the sign functions ensures the symmetric nature 
of the total wave function. In the case of real one-particle eigenfunctions - as is the case 
now - instead of multiplying the determinant with the product of the sign function we 
can take the absolute value of the determinant. Similarly the first excited many body 
state is obtained by replacing the last row in the determinant in Eq. (Tj0) by functions 
of the iV-th excited state of the oscillator. Then, expanding according to the last row we 
get 



1 N 

7^7 E (~l) N+j ~ 1( PN D 3 
'N\j=i 



(6) 



where Dj denotes the (N,j)-ih minor of the determinant in \Ev Due to the orthogonality 
of the single particle functions we obtain 



N-l JV-2 
Qo( x ) = E Wn{x)\ 2 , Q X (x) = J2 l^n^)! 2 + I V^iV (^) | 2 



(7) 



n=0 



n=0 



which can be quickly calculated using a summation formula for orthogonal polynomials 



21] yielding 



-X 2 /I 2 

0o(x) = 2N£ ^ {N _ 1} , (H N ^(x/£)H' N (x/l) - H' N ^(x/£)H N (x/£)) (8) 

and similarly for Q\. The time dependent cross term contains the product of two determi- 
nants both of which can be expanded by their last rows and we obtain 

1 JV N 

*o*i = T^E(-irWi(^ x £(- 1 V -rxlrOP, (9) 
iV! i= i j=1 

Because of the orthogonality of the single particle functions we obtain 

eoi (x) = <PN-i(x)<p N (x) (10) 

which is exact in the one dimensional Tonks model with harmonic trapping. 

Figure Q] shows the time dependent density for 10 atoms in comparison with the ground 



state average density 17| g(x) = g (l — (x/xt) 2 ) 1 ^ 2 , where g = \/2N / (tt£), and xt = 
y2N t is the ground state ID radius of the system. g(x, t) shown here is the fundamental 
swinging mode of the system, and as the number of atoms increases, it approaches the 
average density, while the number of spatial oscillations increases. 

We now consider the signatures of these oscillations in the time-dependent transmission 
of a weak CW field S(x, t), which is sent through the condensate along its axis. This field 
is assumed to be far-detuned from a resonant transition of the atoms of the condensate, 
and will create a polarization density along the sample. The latter can be given as 
V = Nd 2 £/h(/S. — ij), where M is the volume density of the atoms, d is the dipole matrix 
element and 7 is the width of the transition in question, while A = u> — ojl is the detuning 
between the resonant transition frequency ujq and the carrier of the probing field 10 1. We 
write Af = jV a 0(av^y, where M a is the average number of atoms in unit cross section. In 



real experiments % |9j, ll0|] one has a lattice of pencil shaped samples, therefore the actual 
value of M a is the inverse of the cross section of one such "pencil" . (We do not consider 
here the effect of light propagating between these pencils.) The polarization, V leads to 
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FIG. 1: Plot of the time dependent particle density g(x, t) for 10 atoms as function of x measured 
in units of the radius xt- The thin solid line (red in color) is at t = 0, the thick solid line (blue 
in color) is at t = tt/2u and the thin dotted line (orange in color) is at t = tt/uj. We plot the 
average density g(x) for comparison as a black dashed line. 

a space and time dependent susceptibility and index of refraction 



As the time dependence of the harmonic trap is very slow with respect to the frequency 
of the optical fields, instead of the wave equation we shall solve the one dimensional 
amplitude equation for the electric field: 



where k = 2tt/\q is the corresponding wave number in vacuum. Here E(x,t) is the 
temporally slowly varying amplitude of the full electric field: S(x,t) = E(x,t)e~ tuJLt . The 
solution of this equation with a given incoming plane wave from the negative x direction 
shall yield the transmitted wave E tv (t)e tkx for x 3> xt, as well as a reflected wave at 
x <C — Xt- The results of a numerical solution will be discussed below. In order to get 
a better insight into the nature of the problem we also present an approximate analytic 




(11) 




(12) 
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solution to (TT2"]) by a 2nd order WKB approximation obtained for the forward propagating 
wave as 

E,(x ' t) = 7kj) E " exp H nix '' t,dx '-'W-'£ ( ^W) ix ) ' ( 13) 

where E is the incident field amplitude at Xo, far before the condensate. The transmission 
coefficient of the condensate depends only on time for x ^> xt, i.e. far beyond the 
condensate: 



E t {x,t) 



T(t) = ^§^- 2 = exp (-2k £ Im n(x', t)dx' + ±- J\m ^^- dx^j (14) 



since n(x,t) = 1 for x ^> xt- In case of an off-resonant external field we can expand the 
index of refraction given by (fTTl) under the integral up to second order as 

n(x, t) = l + Pq{x, t)/2 - p 2 g 2 {x, t)/8 (15) 

with p = Af a d 2 / (e H(A — ij))- In the following we consider the light source at x = — oo, 
and a photodetector at x = oo. Substituting the second order expansion of the refractive 
index into the formula for the transmission we have: 

/ k Tm B 2 f°° Tm B 2 f°° \ 

Tit) = exp -Nk Jm /3+ P / g 2 (x,t)dx + / (g'(x, t)) 2 dx) (16) 

V. 4 16/c J.oo J 

We write g(x,t) = r (x)+ri(x) cos {ut + a), where r (x) = \co\ 2 Qo(x) + \ci\ 2 Qi(x) is an even 
function and r x (x) = 2|c | |ci|£» i(a ; ) is an odd function of x, and a denotes the relative 
phase of cq and c\. Substituting this into ([TBI) and using the parity of r and ri, we obtain 
the following formula for the transmission: 

T(t) = T N exp [(cos (2(ujt + a))}, (17) 



where 

T N = exp (-Nk Im/3 + ~ Im(/3 2 ) J™ (V 2 (x) + ^r 2 (*)) rfx + ^ f (r* + l(r? ) dx 
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does not depend on time, a = arg(co — cj), and 



C = {2\c,\\c l \) 2 lm{{3 2 ){kR 1 {N)/% + R 2 {N)/Z2k) 



(18) 



with the integrals 



Ri(N) 



Q 2 m (x)dx ) R%{N) 



(g' 01 (x))dx 



depending only on the number of particles in the sample. R\ and R2 can be easily 
calculated up to several hundreds of atoms with a finite sum expression, which can be 



2l| 



derived using known formulas for products of Hermite polynomials 

We can expand the time dependent factor in (fTTl) into a Jacobi form 



22 1, which is 



directly related to the discrete Fourier transform of the time dependent transmission: 



T{t) = T, 



N 



/ (C)+2^/ s (C)cos(2s(u;£ + a)) 



s=l 



(19) 



where I s {() are the modified Bessel functions. This means that the complex Fourier coeffi- 
cients T of the time-dependent transmisson are proportional to modified Bessel functions 
with the same argument, s = 0, 1, 2, ... : 



f{2su) ~ I s {C)e 



-2sia 



(20) 



is rea l) therefore the phase of T(2u) yields the relative phase of the states \l/o and 
a = — \ argT(2u;), in the case of a positive detuning, (ii is odd for real arguments, 
and the sign of ( is the sign of A.) A well known relation for the Bessel functions 22] 
enables us to calculate £ from the transmission spectrum: 



c 



2/i(C) 



2T (2cu)e 



2 1 a 



Io(C)-/2(C) f(0)-\f(4u)\ 



(21) 



Since ( is composed of quantities which characterize the condensate and its interaction 
with the CW field, the knowledge of £ gives information about these quantities. E.g. if we 
know everything in ( except for the coefficients Co and c\, then a measurement of the time- 
dependent transmission yields the value of |co||ci| and the relative phase a which (using 
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FIG. 2: Plot of the electric field amplitude E(x/xT,t = 0), resulting from the numerical solution 
of Eq. (|12[) . with thin solid line (red in color), in arbitrary units, as a probing field passing 
through the condensate of 30 atoms with particle density g(x/xT,t = 0), plotted as a thick 
black line. The approximate analytic solution given by Eq. (|13p differs from the numerical result 
by less then the line thickness in this figure. 

| Co | 2 H- | Ci | 2 = 1) allows us to calculate Co and c\ up to a global phase factor. Alternatively, 
the number of particles in the condensate can be calculated from the Fourier coefficients 
of the time-dependent transmission, if the other quantities in ( are already known. This 
means that one can accurately measure the number of particles in a condensate without 
destroying it. 

We illustrate the use of the modulated transmission for this latter case, by processing 
a simulated transmisson signal which we obtain from the numerical solution of the second 
order amplitude equation ( TT2l) . We assume a sample where 87 Rb atoms are trapped in an 
array of pencil shaped condensates 0, [9], 10] containing 30 atoms in a superposition state 



dl]) with cq = C\ = l/v^2- The laser light for the transmission measurement is assumed to 
be detuned with an angular frequency A = 2tt x 90 MHz from the center of the Di line 



(A =794.978 nm), and we use 7 = 1.80647 x 10 7 1/s and d = 1.4651 x 10" 29 Cm |23fl. Fig. 
[2] shows the density g(x, t = 0) and the electric field amplidude E(x, t = 0) for 30 atoms 
obtained from the solution of Eq. f[T21) . Fig. [3] shows the Fourier amplitudes of the time- 
dependent transmission, assuming a longitudinal trap angular frequency oj = 2n x 100 
Hz. A calculation based on the data shown in this figure and using Eqs. (fT8l) and (T2T1) 
reproduce correctly that the sample contains 30 or 31 atoms. 

In conclusion, we have constructed an exact many body superposition state for the 
Tonks gas, exhibiting a time-dependent particle density in a swinging mode. Generaliza- 
tions for superpositions involving higher excited modes are straightforward. The model 
for the interaction of the Tonks gas with a weak laser beam opens the possibility of mea- 
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FIG. 3: Plot of the magnitude of the Fourier coefficients of the time dependent transmission. 

suring the effect of these density oscillations as a weak but measurable oscillation in the 
transmission signal. The approximate analytic formula obtained for this time-dependent 
transmission allows for the calculation of the quantities which characterize the condensate 
and its interaction with the CW field: the coefficients of the many body superposition 
state, or the number of atoms in the Tonks gas could be measured without destroying the 
sample. 

This work was supported by the Hungarian Scientific Research Fund OTKA under 
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